Fluctuation conductivity from interband pairing in pnictides 
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We derive the effective action for superconducting fluctuations in a four-band model for pnictides, 
discussing the emergence of a single critical mode out of a dominant interband pairing mecha- 
nism. We then apply our model to calculate the paraconductivity in two-dimensional and lay- 
ered three-dimensional systems, and compare our results with recent resistivity measurements in 
SmFeAsOo.sFo.2. 



PACS numbers: 74.40. +k, 74.20.De, 74.25.Fy 

The recent discovery of superconductivity in pnictides 
[H has renewed the interest in high-temperature super- 
conductivity. Pnictides share many similarities with 
cuprate superconductors, e.g., the layered structure, the 
proximity to a magnetic phase 0], the relatively large ra- 
tio between the superconducting (SC) gap and the crit- 
ical temperature T c 0,0], and the small superfluid den- 
sity [2]. However, differently from cuprates, the presence 
of several sheets of the Fermi surface makes the multi- 
band character of superconductivity an unavoidable in- 
gredient of any theoretical model for pnictides. More- 
over, since the calculated electron-phonon coupling can- 
not account for the high values of T c @, 0], it has been 
suggested that the pairing glue is provided by spin fluc- 
tuations exchanged between electrons in different bands 
@, 0, S]- Thus, pnictides are expected to be somehow 
different from other multiband superconductors (e.g., 
MgB 2 ) where the main coupling mechanism is intraband 
i. 

This scenario raises interesting questions regarding the 
appropriate description of SC fluctuations in a multiband 
system dominated by interband pairing. The issue is rel- 
evant, because fluctuating Cooper pairs above T c con- 
tribute to several observable quantities, such as, e.g., the 
enhancement of dc conductivity (paraconductivity) and 
of the diamagnetic response as T c is approached jlfjj] . The 
nature of SC fluctuations depends on whether the system 
is weakly or strongly coupled (and on whether preformed 
pairs are present or not) and a wealth of physical infor- 
mation can be obtained from paraconductivity and dia- 
magnetic response, provided a theoretical background is 
established to extract them. 

In this paper, after introducing a four-band model, as 
appropriate for pnictides, we discuss the subtleties re- 
lated to the description of SC fluctuations in a system 
with dominant interband pairing. We then apply our 
results to compute the paraconductivity associated with 
SC fluctuations above T c . We show that when interband 
pairing dominates, despite the presence of four bands, 
there are only two independent fluctuating modes. Only 
one of them is critical and yields a diverging Aslamazov- 



Larkin (AL) contribution to paraconductivity as T c is 
approached, similarly to the case of dominant intraband 
pairing [lj]]. The temperature dependence of the AL 
paraconductivity is the same derived for ordinary single- 



band superconductors [10|, LL2|]. Remarkably, within a 



BCS approach, we recover the AL numerical prefactor, 
which is a universal coefficient in two dimensions, and 
depends instead on the coherence length perpendicular 
to the planes in the three-dimensional layered case. We 
also find that subleading terms with respect to the lead- 
ing AL contribution could distinguish between dominant 
inter- and intra-band pairing. Within this theoretical 
background, we analyze recent resistivity data in Sm- 
FeAsOF 0] and discuss the results. 

At present, ARPES measurements in pnictides 0] con- 
firmed the Fermi-surface topology predicted by LDA. It 
consists of two hole-like pockets centered around the T 
point (labeled a and /?, following Ref. [a]), and two 
electron-like (7) pockets centered around the M points 
of the folded Brillouin zone of the FeAs planes. Moti- 
vated by the magnetic character of the undoped parent 
compound and by the approximate nesting of the hole 
and the electron pockets with respect to the magnetic 
ordering wavevector, we assume that pairing mediated 
by spin fluctuations is effective only between hole and 
electron bands @. Since the (3 band has a Fermi sur- 
face substantially larger than the a band, the {3 pocket 
is expected to be less nested to the 7 pocket, so that the 
interband (3 — 7 coupling is smaller than the a — 7 cou- 
pling , i.e., rj = A/A < 1. The two electron pockets have 
comparable sizes and for simplicity we assume that the 7 
bands are degenerate. Therefore, the BCS Hamiltonian 
of our four-band model is f 131 



q 



is the band Hamiltonian, 



where flg = £ k &,kc] jktr c i>k(T 

c i ktr annihilates (creates) a fermion in the i — a, (3, 71, 72 
band (with the twofold degenerate 7 bands labeled as 71 
and 72), £i,k is the dispersion with respect to the chemical 



2 



potential, <J>i !q = X)k c i,k+qT c i,ki ls the pairing operator 



in the z-th band, and $ n 



+ 3> 7 , n • Since we as- 



sumed that pairing acts between hole and electron bands 
only, we can express the pairing term in Eq. Q] by means 
of $1 = $ 7 and <J> 2 = <I> Q + rj^p. Thus 



= -A 



q 



), (2) 



where $± = (idi$i ± u^^Vv^, wi,W2 are arbitrary 
numbers satisfying wii«2 = 1, and for definiteness we 
take A > 0, which is the case for a spin-mediated pair- 
ing interaction. From Eq. one immediately sees that 
when interband pairing dominates, the interaction is a 
mixture of attraction (for <&_) and repulsion (for $ + ). As 
a consequence, when we perform the standard Hubbard- 
Stratonovich (HS) decoupling of the quartic interaction 
term ^ by means of the HS field <j), 



±A$t$ 



V(j>e 



-\c/>\ 2 /A+^/±T(&<p+h.c.) 



(3) 



the HS transformation associated with the repulsive part 
contains the imaginary unit. As we shall see below, this 
would require an imaginary value of the 4>+ HS field at 
the saddle point, that can be avoided by the rotation of 
the pairing operators $± defined above, via a suitable 
choice of the w\ t 2 coefficients. 

In analogy with the single-band case the effective 
action for the SC fluctuations reads 



s = EEE&* 
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0*(g)$_(g)-^-(g)*L(g) 



where k = (k, is n ), q = (q,iw m ), and e n , oj m are the 
Matsubara fermion and boson frequencies, respectively. 
Integrating out the fermions we obtain the standard con- 
tribution to the action, — i Tr 1°S ^\k' > wnere the trace 
acts on momenta, frequencies, and spins. The elements 



of the matrices A\. k , 


are 
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_(fc - fc') + #+(fc - k')} , 



for 7 = 71, 72, and the [Aj^^i elements contain the com- 
plex conjugates of the HS fields evaluated at (k' — k). 



The q — k — k' — values of the HS fields yield the 
SC gaps in the various bands, within the saddle-point 
approximation 4>±{q = 0) = (f>±. However, due to the 
presence of the imaginary unit i associated with the HS 
field <j> +j in general [^ fc ]2i [-^jad 12 • To recover a 
Hermitian A matrix at the saddle point, the integration 
contour of the functional integral must be deformed to- 
ward the imaginary axis of the 4>+ field. This can be 
avoided if one chooses w\ and u>2 in the definition of 
<&± in such a way that (f> + — 0. In our case, the choice 
4>+ = gives A Q = — zi>2</>_/\/2, A 7 = w\4>-/\/2 and 
A^ = r]A a . Hence, the ratio of the gaps in the two hole 
bands is solely determined by the ratio of the couplings, 
A^/Aq = i). Using w\W2 — 1 we have <P = — 2A Q A 7 , 
and the equations for A Q and A 7 read [131] 



7// A 7 , 



A Q = -A(2n 
A 7 = —A (n a A Q 



(4) 
(5) 



where n, = N, ^° ^[tanh(£; i /2T)]/ J B. i yields the q = 
value of particle-particle bubble when T > T c , Ni is den- 
sity of states of the i-th band at the Fermi level, ujq is the 
cut-off for the pairing interaction, and Ei = ^/^ 2 + Af. 
Hence, our four-band model for pnictides effectively re- 
duces to a two-band model, with one electron-like and 
one hole-like effective band. Indeed, defining Ai = A 7 , 

a 2 = A Q , and n x = 2n 7 , n 2 = n Q + ^itj, E qs . Q-© 

recover the standard two-band expression Ai = — An 2 A 2 
and A2 = — AniAi. Since A1/A2 = —w\/w2, it also fol- 
lows that w\/wl = U 2 /Ui at T < T c . 

Let us now discuss the SC fluctuations for T > T c , 
where Ai, 2 = 0. To derive the equivalent of the standard 
Gaussian Ginzburg-Landau functional [ltj, we expand 
the action up to terms quadratic in the HS fields, Sg = 
E t « E g 4>t{9)L^{q)(t) K {q), where i,k = ± and 



A" 



1 -^eff{q) -iSeff(q) 
-iEeffiq) A-i+ILeffiq) 



(6) 



with 



n e/ /(?) 



1 , 



ln 2 (g)], 
n 2 (<?)]. 



s e//(<?) = -[w-Jli{q) - w 2 



The critical temperature is determined by the condi- 
tion detL~ 1 (q = 0) = 0. In the BCS approximation 
ni(0) = ATjln(1.13 u>q/T c ) and, in agreement with the 
Eqs. Q-©, we obtain N e// Alog(1.13 uj /T c ) = 1. Here, 
the parameter N eff = y / N 1 N 2 = y/2N^(N a + r) 2 Np) 
plays the role of an effective density of states. To com- 
pute the fluctuation contribution to the various physical 
quantities at T > T c , we evaluate Eq. © at leading order 
in q (hydrodynamic approximation), using the standard 
expansion of the particle-particle bubble for a layered 
system, n,(g) w H^O) 



c i,\\1\\ 



with 



3 



Q\\ = + 1y an d the z axis perpendicular to the layers. 
In the BCS approximation, e.g., ji = %Ni/ (8T). We omit 
the lengthier but standard BCS expressions for c^y and 
Ci iZ (see, e.g., Ref. [13]), that will not be explicitly used in 
the following. By making for T > T c the same choice pre- 
viously adopted for T < T c , wfTLi(0) = 112(0), we sim- 
plify the structure of the fluctuating modes. Indeed, with 
this choice, the off-diagonal terms in Eq. ([6]) yield con- 
tributions beyond hydrodynamics and can be neglected. 
Therefore, the leading SC fluctuations are described by 



ni- + v(q) 

m+ — v{q) 



(7) 



where m± — A -1 ± v / lLTl2 are the masses of the collec- 
tive modes, v{q) = c\\qj + c z q 2 z + j\u m \, cy = {w\c lf \\ + 
w%c 2 ^\)/2 is the stiffness along the layers, c z — (w 2 Ci jZ + 
w|c2 lZ )/2 is the stiffness in the direction perpendicular 
to the layers, and 7 = (w^ji + WYi%)j1 is the damping 
coefficient. In the BCS case m± — N e f /[ln(1.13 ujq/T c ) ± 
ln(1.13w /T)] and 7 = nN ef f/(8T). 

Having deduced the hydrodynamic action of the SC 
fluctuations, we can calculate the paraconductivity along 
the lines of Ref. By inspection of Eq. (7j), one can 

see that only the <fi- mode becomes critical at T c [i.e., 
to_(T c ) = 0, m + (T c ) = 2A _1 ], thus giving a diverg- 
ing fluctuation contribution to various physical quantities 
when T — > T c . The leading contribution to the current- 
current response function along the layers is 



S X (n e ) = 4e 2 T^ C fg 2 L__(q,a; m ) J L__(q,w„ 



whence the paraconductivity 8a al — [Irn£x(f2)/^]n_>o 
is obtained, after analytical continuation of the Matsub- 
ara frequency fi^ to the real frequency Q. Therefore, 
the same expressions known for a single-band supercon- 
ductor are found, although with the effective parameters 
m_, c z , 7. The paraconductivity along the layers is inde- 
pendent of the in-plane stiffness cy , as guaranteed by the 
same gauge-invariance arguments discussed in Ref. [HI 
for a single-band superconductor. The leading contribu- 
tions to paraconductivity along the layers in three and 
two dimensions (3D and 2D, respectively) take the AL 
form El 



1 



32/16 y/l 

e 2 1 



(8) 
(9) 



where £ z — yCz/l ls the correlation length in the di- 
rection perpendicular to the layers, d is the distance be- 
tween layers, and e(T) = nni-(T)/(8jT c ) is the dimen- 
sionless mass of the critical collective mode. We point out 
that the above expressions are general within a hydrody- 
namic description of the collective modes and do not rely 



on any particular assumption about the pairing strength 
When the BCS expression for the IT bubbles holds, 
m_(T) = N e ffln(T/T c ) and the dimensionless mass ap- 
pearing in Eqs. |(8])-((9]) is simply e = \og(T/T c ). 

The calculated AL paraconductivity expressions ([8])- 
§2§ can be now compared with the existing results for 
the pnictides. As it has been widely discussed in the 
context of cuprates jl6], for weakly-coupled layered ma- 
terials the SC fluctuations usually display a 2D-3D di- 
mensional crossover as T c is approached. However, the 
interlayer coupling has a different relevance in the various 
families of cuprates, with substantial 3D behavior [Eq. 
JH])] in YBa2Cu306+a; samples, while more anisotropic 
Bi 2 Sr 2 CaCu208 or La^^Sr^CuC^ compounds show 2D 
fluctuations [Eq. the 2D-3D crossover being too 

close to T c to be clearly observed. Such a systematic sur- 
vey has not yet been performed in the case of pnictides, 
due also to the limited availability of clean single crys- 
tal. Indeed, the analysis of the 2D-3D crossover might 
be biased in polycrystals by the distribution of critical 
temperatures and by the mixing of the planar and per- 
pendicular directions. 

Having in mind such limitations, we attempt the anal- 
ysis of paraconductivity in a SmFeAsOo.sFo.2 sample 
[4j with T c fa 52 K. To determine the contribution of 
SC fluctuations to the normal-state conductivity, 5a = 
p~ x — p^y 1 , we need to extract the normal state resistivity 
Pn from the data. Owing to the diverging conductivity, 
the precise determination of the normal-state contribu- 
tion is immaterial near T c , but becomes relevant for larger 
values of e. We fitted the resistivity at high temperature 
(in the range between 279 K, the highest available tem- 
perature, and about 200 K) checking that the qualitative 
results were stable upon small variations of this range. 
We used a quadratic fit pjv — a + b(T — To) + c(T — To) 2 , 
with To = 279 K, and found that the resulting paracon- 
ductivity is roughly two orders of magnitude smaller than 
the 2D AL result [Eq. ©]. The slope in a log- log plot is 
in agreement with the 3D power law [Eq. J5])]. The plot 
in Fig. Q] clearly shows that the 3D behavior extends over 
two decades of e. The fitting parameters a = 1700 p,Q cm, 
b = 21.9/zficm/K, and c = 0.018^ cm/K 2 were used. 

The fitting with Eq. ([8]) allows us to determine the 
precise value of T c and the coherence length £ 2 of the SC 
fluctuations in the direction perpendicular to the FeAs 
layers. We find T c = 51.4 K and £ z = 19 A. One can 
also see that substantial SC fluctuations persist up to 
e w (T - T c )/T c w 0.4, i.e., up to 20 K above T c . This 
fluctuating regime is therefore much larger than in con- 
ventional 3D superconductors and, despite the 3D behav- 
ior of the paraconductivity, calls for a relevant character 
of the layered structure and for a small planar coher- 
ence length. We point out that at even larger values of 
e the paraconductivity drastically drops, in analogy with 
what found in cuprates (see, e.g., Ref. [16], and refer- 
ences therein). The way paraconductivity deviates from 
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0.1 




FIG. 1: (Color online) Comparison between the experimental 
paraconductivity for a SmFeAsOo.8Fo.2 sample studied in Ref. 
[J] (black circles) and the 2D (dashed line) and 3D (solid line) 
expressions of the AL paraconductivity [Eqs. I|8I9II ]. For the 
3D paraconductivity a coherence length £ z = 19 A has been 
used, while for the 2D case the structural distance between 
layers d = 8.4 A has been inserted. 

the AL behavior in multiband systems also depends on 
the role of the other (non-critical) collective modes. In 
particular, it can be shown [ijj that, when the intraband 
pairing is equally dominant in all bands, the paraconduc- 
tivity mediated by the non critical modes may become 
sizable, and the experimental data should approach the 
pure AL contribution of the critical mode from above. 
This is not the case when the dominant pairing is in- 
terband and therefore it is not surprising that the data 
for the pnictide sample analyzed in this paper always lay 
below the AL straight line in Fig. 1. 

In conclusion, we investigated the occurrence of SC 
fluctuations in a multiband system where interband pair- 
ing dominates, as appropriate for pnictides. In contrast 
to the case of dominant intraband mechanism (as, e.g., 
in MgB 2 [LH), in the present situation the HS decou- 
pling must be accompanied with a proper rotation of the 
fermion fields which guarantees a Hermitian saddle-point 
action below T c . The same rotation leads to a straight- 
forward decoupling of the Gaussian fluctuations above T c 
in the hydrodynamic limit. Thus, despite the apparent 
complexity of the multiband structure in pnictides, we 
demonstrate that the AL expressions for paraconductiv- 
ity stay valid not only as far as the functional temper- 
ature dependence is concerned, but also regarding the 
numerical prefactors. While in the BCS 2D case the 
prefactor stays universal, in the 3D case the only dif- 
ference is that a suitable redefinition of the transverse 
coherence length has to be introduced. With this equip- 
ment, we considered the experimental resistivity data of 
the SmFeAsOo.sFo.2 sample studied in Ref. [4J|, finding 
that here fluctuations have a 3D character and extend far 



above T c . Recently, an experimental work[18] on fluctu- 
ation conductivity in pnictides confirmed the wide fluc- 
tuating regime, even though fluctuations seem to have 
2D character. Thus, further experiments are in order to 
confirm the nature of fluctuations in pnictides, and to 
assess the relevance of Cooper-pair fluctuations in these 
new superconductors. 
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